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We provide evidence for the mapping of critical spin-1 chains, in particular the SU(3) symmetric bilinear- 
biquadratic model with additional interactions, to free boson theories using exact diagonalization and the density 
matrix renormalization group algorithm. Using the correspondence with a conformal field theory with central 
charge c = 2, we determine the analytic formulae for the scaling dimensions in terms of four Tomonaga- 
Luttinger liquid parameters. By matching the lowest scaling dimensions, we numerically calculate these field- 
theoretic parameters and track their evolution as a function of the parameters of the lattice model. 


I. INTRODUCTION 


renormalization group [23] methods [15, 24-27]. 


Given a strongly correlated quantum system, an important 
step towards understanding it is to determine its basic prop¬ 
erties, such as the presence or absence of a gap, the presence 
of spontaneous symmetry breaking, etc. One then asks for 
more specific information, and ultimately, the complete de¬ 
scription of the underlying low-energy physics. Quite often, 
this characterization involves determining an effective field 
theory. Examples include topological field theories describ¬ 
ing the full braiding statistics in gapped quantum systems, and 
conformal field theories (CFT), describing the set of indepen¬ 
dent critical exponents in gapless systems. Obtaining these 
conformal exponents is important because close to the critical 
point, the power law behavior of physical quantities like mag¬ 
netic susceptibility is governed by them. These dimensions 
complement the knowledge of the central charge, denoted by 
c, in determining the universal long-distance behavior of the 
theory. 

In recent times, several probes, such as the entanglement 
entropy (EE), Renyi entropies and entanglement spectrum [1- 
9], have been devised to explore the above mentioned proper¬ 
ties. A central component of all these measures is the ground 
state reduced density matrix, calculated for a finite region of 
space and obtained by tracing the full density matrix over the 
other degrees of freedom. For example, the finite-size scaling 
of the EE in one-dimensional critical systems provides a pre¬ 
cise estimate of the central charge of the corresponding CFT. 
More sophisticated ways of using reduced density matrices 
also reveal information about the low-energy scaling dimen¬ 
sions and operators [10-15]. 

For one dimensional (ID) critical systems, the theory of 
Tomonaga-Luttinger liquids (TLLs) [16-19] has been remark¬ 
ably successful at characterizing their low-energy physics. 
There has been additional validation on the experimental 
front, at least qualitatively; several realizations, ranging from 
carbon nanotubes [20, 21] to semiconductor wires [22], of 
TLL physics have been found. Quantitative estimates of 
the scaling dimensions, velocity, and Luttinger parameter for 
model Hamiltonians have been made with analytic solutions 
or numerically, with exact diagonalization and density matrix 



Figure 1. A representative configuration of hard-core bosons on a 
three leg ladder with periodic boundary conditions in the transverse 
or rung direction, with Hamiltonian given by Eq. (la). The hopping 
in the length and transverse {ty) directions respectively and the 
correlated hop {K) on the square plaquette have been indicated by 
arrows. For = 0 and at 1/3 filling, the low energy model involves 
configurations with only with exactly one boson per rung. The three 
configurations per rung, one of which has been enclosed in a dotted 
rectangle, when appropriately Fourier transformed, are mapped to a 
spin-1 basis. The resultant spin-1 Hamiltonian has the form Eq. (2a) 

Most theoretical studies have focused on the single compo¬ 
nent TLL, which directly corresponds to a c = 1 CFT, and 
which now appears to be a fairly well understood case [15, 
18]. In contrast, there are few general results for the c = 
2 case, despite the existence of systems with this prop¬ 
erty [28, 29]. This is partly attributed to the TLL theory for 
c = 1 being completely described by a single dimensionless 
parameter, whereas the c = 2 theory requires four dimension¬ 
less parameters. An important open question is that there is 
no established method to extract TLL parameters for a given 
lattice model. Given the history of the TLL, it appears to us 
that this situation is quite unsatisfactory and incomplete. 

In special cases, a c = 2 CFT can be understood as a tensor 
product of two c = 1 CFTs; for example, the ID Hubbard 
model has two TLL parameters, one for spin and the other for 
charge. In this paper, however, we will demonstrate a TLL 
parameter extraction procedure for a c = 2 CFT where such 
a decomposition does not apply. Several conceptual and prac¬ 
tical questions arise here; including which measures must be 
calculated to estimate them and whether they are numerically 
accurate enough to validate or refute a given field theory. Our 
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paper addresses these questions and highlights an interesting 
application of relatively new ground state entanglement based 
metrics, such as the mutual information. However, before 
considering a specific problem to demonstrate our ideas for 
c = 2 CFTs, we mention physical examples where this situa¬ 
tion occurs. 

One way to realize a multi-component TLL is to couple sev¬ 
eral TLLs [30, 31]. The most natural geometry for doing this 
is a ladder (or tube), a quasi-one dimensional system made up 
of one dimensional legs which are additionally coupled in the 
transverse or rung direction, with open (or periodic) bound¬ 
ary conditions. Fig. I shows an example with three legs, rel¬ 
evant for modelling quasi one-dimensional compounds such 
as [(CuCl 2 tachH) 3 Cl]Cl 2 [32] and CsCrF 4 [33], and to which 
recent theoretical works [29, 34-36] have been devoted. 

Following the work of Ref. [34] and as is schematically de¬ 
picted in Fig. (I), our starting point is a system of hard-core 
bosons on a three leg tube, governed by the Hamiltonian, 

H = TJhop + Hk (la) 

Ffhop — fa; ^ ^ h.c. ty ^ ^ b^b^-^-y h.c.(lb) 

r r 

Hr = blhx+xblj^^j^^hr+y + h.c. (Ic) 

r 

where ty and are the hoppings along the transverse {y) and 
length (a:) directions respectively, iF is a correlated exchange 
on a square plaquette. The phase diagram of this model is 
expected to be quite rich; here we only consider the case 
of tx,ty —i' 0 with 1/3 filling of bosons. In this parameter 
regime, the low-energy theory of this model involves only one 
boson per rung (column) allowing three distinct states on it; 
the number per rung can not change because of the absence of 
hopping in the x direction. 

On Fourier transforming bosonic creation operators along 
the y direction, a new basis at every x location is 
defined as, |0)^ = ^ ^ac) and 

(^1,0where w = 
exp (i27r/3). The three components can be thought of as those 
corresponding to a pseudospin-1 object, leading to the effec¬ 
tive spin-Hamiltonian of the form [34], 

iF = S* • S, + ^ (S, • ) + Q, + Q^a) 

(*/> (*i> 

Qx{y) = Qx(y) ^ 
i 

where is a spin-1 operator living on site i, while and qy 
are scalars. K will be set to 1 throughout and thus all energy 
scales in this paper are in terms of this unit. Ux and Uy are 
3x3 matrices in the Sz basis (ordered as —1, 0,1), and are 
given by, 

/w-i00\ /001\ 

Ux=[ 0 10, [/j, = 1 0 0 , (3) 

yoow/ \oio/ 


The Qy term is non-zero with qy = —1/3, but there is no 
Qx term i.e. qx = 0. Physically, the term Uy’^U^ models a 
correlated cyclic permutation of neighboring spins. However, 
diagonalizing Uy, i.e. performing a similarity transformation 
by the matrix. 




(4) 


preserves the combined SU(3) symmetry of the first two 
terms in (2a) and converts the Qy term into the Qx term be¬ 
cause SUyS"^ = Ux- For presentational purposes, we have 
shown both terms in Eq. (2a); this generalized model has been 
previously introduced in the literature as the quantum torus 
chain [37]. 

For qx = Qy = 0, this model is the analytically solvable 
Lai-Sutherland model [38-40], which serves as a useful guide 
for checking our calculations. Since Qx and Qy are related 
by a SU(3) unitary transformation; studying the model with 
qx non zero and gj, = 0 is equivalent to the case with qx = 
0 and qy non-zero. We set qy = 0 throughout this paper, 
and leave the more general case for later exploration. Finally, 
we note that a generalized version is the bilinear-biquadratic 
model [41], whose phase-diagram includes a gapless phase 
and the gapped Haldane phase [42, 43] and which has been 
experimentally realized in LiVGe206 [44]. 

We now discuss the organization of the remainder of the 
paper. In Sec. II, we discuss how the low-energy theory of the 
spin-1 model (2a), motivated above, is mapped to a field the¬ 
ory using bosonization techniques. We then develop the an¬ 
alytic formulas for the scaling dimensions of the low-energy 
theory in terms of the TLL parameters: these formulae are 
generalizations of those known for the c = 1 case [45]. Lor the 
particular case of parameters of the spin-1 Hamiltonian (2a) 
{qx > 0, qy = 0), these formulae show the explicit depen¬ 
dence of the TLL parameters on the microscopic model pa¬ 
rameter. In Sec. Ill, we provide numerical evidence for the 
connection between the low energy theory of the spin chain 
and the CLT for c = 2, by calculating the lowest two scaling 
dimensions with exact diagonalization (ED) and the density 
matrix renormalization group (DMRG). The TLL parameters 
obtained are tracked as a function of the microscopic param¬ 
eter qx- Linally in Sec. IV, we conclude by discussing gen¬ 
eralizations of our method and the prospective applications to 
other systems. 


II. MAPPING SPIN-I LATTICE MODELS TO EREE 
BOSON THEORY 

A. Symmetries 

In this section, we develop a continuum field theory de¬ 
scription for the lattice Hamiltonian (2a), by closely follow¬ 
ing Refs. [40, 46], wherein more details are spelled out. Lor 
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a start, symmetry properties of the Hamiltonian (2a) are de¬ 
scribed here piece by piece. To this end, the spin-1 part of the 
Hamiltonian (2a) can be written (up to a constant factor) in a 
manifestly SU(3) symmetric way in terms of 3 x 3 elemen¬ 
tary matrices with one on row a and column /3 and zero 
everywhere else as 

^SU(3) = (5) 

{ij) a,/3=0,l,2 

For convenience, this Hamiltonian can be represented in 
terms of fermionic operators using = c^^cpi, with 

the constraint Ea=o,i.2 = J 2 a=o.i, 2 ^ai = 1 at 

each site i. The constraint ensures that the operators have 
the same commutation (anticommutation) relations and act 
on Hilbert spaces of the same dimensions. The Hamilto¬ 
nian iFsu( 3 ) conserves the particle numbers Ni — Nq and 
N 2 — Nq, where Na = Defining the dual basis by 

c„ = 3“^/^ Sa=o for '' 2 - = 0,1, 2, and the correspond¬ 
ing particle numbers as 1V„ = the Hamiltonian 

iFsu( 3 ) ^Iso conserves the dual particle numbers Ni — Nq and 
N2 - Nq. 

On the other hand, the perturbation in (2a) can be writ¬ 
ten as 

= ( 6 ) 

i,ct 

The Hamiltonian Tfsu( 3 ) +Ox conserves the particle numbers 
Ni — Nq and N 2 — Nq, and it conserves the dual particle 
numbers Ni — Nq and N 2 — Nq (mod 3). 


shown to be marginally irrelevant in the RG sense. While this 
term must be retained to evaluate quantitative finite size loga¬ 
rithmic corrections, here we simplify the analysis by working 
directly in the conformal limit. Instead, the finite size cor¬ 
rections will be reintroduced only at a later stage, when com¬ 
paring the analytic results with numerical calculations. Thus, 
with this simplification, the critical theory is, 

Hsvo) « i?wzw = TTI dxY,[J^p4^ + JlpjLl 

OL,j3 

( 10 ) 

We note that this is a C/(3) Wess-Zumino-Witten (WZW) 
model and thus contains an SU(3)j^ and (7(1) part [47]. The 
(7(1) piece is precisely the charged mode which is gapped 
and will be dropped later. (N.B. the above procedure is better 
described and applied, instead of dealing with the SU(3) Lai- 
Sutherland model, by starting with the Hubbard type model 
H = -f E(ii)a[claSa + ^ withoUt 

constraint Uai = 1. This constraint is in fact generated 
dynamically and this model reduces to the SU(3) symmetric 
spin model when expanded in t/U) 

Applying the same reasoning as above one deduces the con¬ 
tinuum approximation 

dxY,[{JRS + {Jlc.f]. ( 11 ) 

a 

where again we have dropped the terms that only depend on 
the charged mode. 


B. Continuum theory 


C. Abelian Bosonization 


The low-energy effective field theory for the Hamiltonian 
i7su(3) can be developed by noting that at low energies, only 
excitations close to the Fermi points kp = 7r/3ao (where uq 
is the lattice constant) propagate. Thus we can approximate, 

Caz « ^La(Xz)] ■ (7) 

Substituting this in the Hamiltonian and dropping oscillatory 
terms, the low energy theory can be written in terms of the 
(7(3) currents. 


TO 


'4’rIl'^r,lp, 


( 8 ) 


as 

Hsvo) = 7rvp I dxY,[Jlp4^ + Jlp4^+2J’^p4^ 

- 2 cos(2fcFao) p], (9) 

where uf is the Fermi velocity, which will be set to 1 hence¬ 
forth. The last term depends only on the charged modes, 
and which are gapped, while the second term can be 


Introducing holomorphic and antiholomorphic coordinates 
for 1 h- 1 d space-time z = —i{x — t), and z = i{x + t), the 
time evolution of the fields factorize nicely so that the fields 
with an R{L) subscript depend only on z{z) respectively. The 
continuum Fermi fields in Eq. (7) can be bosonized as follows 


i’pL — 


27rao 


„—iV4^<Pl3L 


= 


27rao 




( 12 ) 


where (pp represents the holomorphic part of a free boson 
field. We focus on the holomorphic parts of the theory (drop¬ 
ping the L subscript) with similar formulae for left moving 
fermions in terms of the anti-holomorphic part of the free bo¬ 
son field. We have introduced normal ordering of an operator 
O, denoted by : O :, which must be used when two fields at the 
same point are multiplied together. Usually when bosonizing 
more that one species of fermions, one introduces Klein fac¬ 
tors to ensure that different Fermi fields anticommute. These 
Klein factors have been ignored here since they are not dy¬ 
namical and do not play a role in the Hamiltonian which is 
mainly what we are interested in here. 
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A key identity, which can be regarded as the inverse of Eq. 
( 12 ), is 

= -^d(t)a{z), (13) 

VTT 

where d denotes a derivative with respect to z. To understand 
the splitinto SU(3) and U(l) WZW theories mentioned above 
we introduce the SU(3) and U(l) currents given by 

r = Y, j = Y. (14) 

q :,/3 a 

where r“ are generators of the SU(3) algebra. The C/(l) 
piece in the boson language satisfies 

J = —^ + d(j)i + d(j) 2 ) ■ (15) 

y/TT 

The SU(3) currents associated to the Cartan sub-algebra are 

cx cx d(j)Q — d4>i, 

cx ~ 2 ' 0 ^ 1'02 oc 84)0 + d4>i — 2902- 

(16) 

So we can make an operator product expansion (OPE) pre¬ 
serving orthogonal change of basis to introduce 4 >o^i .2 as 

00 = (00 + 01 + 02)/v^, (17) 

01 = (00 - 0l)/V^, 

02 = (00 + 01 - 202)/v/6. 

In this basis the dynamics of the charged mode is now encoded 
in the single boson field 0o. Therefore dropping the charged 
mode corresponds to setting 0o = 0. This is indicated with an 
arrow in the equations below. In this basis some of the SU(3) 
currents associated with the Cartan subalgebra are simply (up 
to a constant factor) 90i, 902, while those associated with a 
choice of simple roots for SU(3) are 

cx 0j0i CX ^ 

J“" cx 0^02 oc ^ 

Qci and a 2 together with a third root a 3 are given by, 

o.=(l, 0 ),aj= ( 1 ,-/?). (18) 

All other roots of SU(3) can be obtained as integer linear com¬ 
binations of 0.1 and 02 , for example 03 = Oi — 02 - Sim¬ 
ilarly, the vertex operators associated with all other roots can 
be obtained from operator products of and 9“^. This con¬ 
struction gives precisely the vertex operators obtained in the 
purely bosonic construction of SU(3)j^ where the boson fields 
01 and 02 are compactified on the root lattice of the SU(3) 
algebra [48]. All proportionality constants can be fixed by a 
choice of normalization of the SU(3) generators. 


We now obtain the purely bosonic description of the gapless 
degrees of freedom of the spin-1 chain. The key identities are 

Y.j:j: = -Y.d4>ad4>^ 

a Oi 

-(90i90i -I- 902902), (19) 

E = - E 

-)■ -(90i90i -I- 842842 )- (20) 

Using these, we obtain the main results of this section 

77 su( 3 ) ^ J dx{ 84 i 84 i + 902002 + antilhol), 

Qx ^ J da;(90i90i-I -902002 + antilhol). ( 21 ) 

where antihol denotes the antiholomorphic part. Note that the 
tildes have now been dropped; the effective Hamiltonian of 
the spin-1 chain is now written in terms of SU(3) 1 boson fields 
01 and 02 . We note that these quantities are all non-negative 
since we have for any field 0 , 

9090 -f 9090 = - ^ {8x49x4 + 8t4dt4) ■ (22) 

D. General c = 2 Boson theories 

In the previous section, we derived the low-energy effec¬ 
tive Hamiltonian that should capture the critical dynamics of 
the spin-1 chain with the Qx perturbation. The low-energy 
effective theory consists of two compactified boson fields and 
has the central charge c = 2. To put the effective theory in a 
general context, we discuss in this subsection a generic two- 
component boson theory with c = 2 . 

Eor the case of the single-component TEL, the landscape of 
the theory (often called “moduli space”) is well understood. It 
is characterized solely by a single parameter, the Luttinger pa¬ 
rameter K or the compactification radius R of the boson field. 
There is a boson-vortex duality in (iH-l)d (also known as “T- 
duality”) which relates the two regions K > \ and K < 1. 
These regions are separated by the self-dual point K — \ 
where SU(2) symmetry is realized. With orbifolding, theory 
space for c = 1 is described in terms of two axis, each de¬ 
scribing the ordinary free boson theory (the single-component 
TEL) and its orbifolded counterpart, together with a few “ex¬ 
ceptional cases” [49]. 

On the other hand, the moduli space for the c = 2 theories 
is more complicated. Eor a start, let us consider the action in 
1 h-1 d space-time for two bosonic fields 

S=^J dxdt iGab8^X‘^8^X^ + Babe^,8^X‘^8,X^), 

( 23 ) 
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where fj,,v = 0,1, and G and B are a symmetric (non degen¬ 
erate) and antisymmetric 2 by 2 real matrix, respectively. 

-h 27r. (24) 

The corresponding Hamiltonian is given by, 

H = f dx + antihol), (25) 


torus of modulus t = ti + iT 2 i.e we compactify Euclidean 
space-time as a: ~ a; -f 27r and (a:, t) ~ {x, t) + 2 'xt 2 ). 

One can quantize using path integrals, the path integral yields 
a sum over instanton sectors. We can write in each instanton 
sector 


where 


vCL _ 

^n,w '^n,iu,cl 




(29) 


Observe that the parameter B does not enter into the Hamilto¬ 
nian: it is a topological term. However, it affects the canonical 
commutation relations and hence the spectrum. 

There are thus four independent parameters, Gn, G 12 , G 22 
and Bi 2 , characterizing the c = 2 action (23), as opposed to 
the c = 1 TLL parameterized by a single parameter. (As in 
the case of c = 1, one can consider various orbifolds of the 
two-component boson theory (23), leading to an even richer 
moduli space or phase diagram [50]. ) 

For the case of c = 1 TLL, the duality relates the large 
and small compactification radius (the Luttinger parameter). 
Similarly, there is a group of duality transformations acting 
on the four parameters, and different values of G and B do 
not necessarily correspond to different spectra [51, 52]. To 
unveil this duality group, it is convenient to trade the four real 
parameters in G and B for two complex parameters ^ and p 
as follows 


f _ Gi2 
P = Bi 2 - 


, v/detG 

* ’ 
(J22 

- iv/det G. 


(26) 


These two parameters can be acted upon by independent 
SL{2, Z) transformations which for ^ is given by 




-I- 6 

+ d’ 


(27) 


where a, b,c,d G Z, ad — be = 1. There is a similar inde¬ 
pendent transformation for p. These transformations change 
the parameters G and B but lead to the same spectrum. Effec¬ 
tively the target space of the boson fields corresponds to two 
tori, which are left invariant by SL{2, Z) x SL{2, Z) transfor¬ 
mations. There are two further discrete transformations that 
leave the spectrum invariant: 


(C) p) (Pj ^)) (^) p) (~P) ~C)- (28) 


When Bi 2 = G 12 = 0, we have a product of two c = 1 theo¬ 
ries. In this case, the first transformation sends G —G“^, 
which corresponds to two independent duality transforma¬ 
tions for each c = 1 theory. Fig. 2 depicts a portion of the 
space of theories in the ^ = p plane together with some points 
of enhanced symmetry. We anticipate that these theories cap¬ 
ture the critical behavior of the gapless degrees of freedom of 
spin-1 chains such as the model in Eq. (2a). 

To deduce the spectrum of these bosonic theories we switch 
to Euclidean signature t —)■ —it. We take space-time to be a 


„ , . „ {n°‘— w°'Ti)t 

+ (30) 

T 2 

is a classical solution that winds n and w times along the two 
non trivial cycles on the torus. The partition function is 

^ j (31) 

n,w 


where is the classical action evaluated on shell for Eq. 
(30), and the quantum path integral is over a continuous un- 
compactified variable Xg. The second term in Eq. (23) is 
a total derivative for periodic functions Xq and can be ne¬ 
glected. Thus the integral over Xq yields the determinant of 
the quadratic differential operator appearing in Eq. (23) which 
is just the Laplacian in the spacetime index times the G matrix 
in the internal index a. After applying a Poisson resummation 
in n for the classical contribution one finds 


Z = 


( V 

Vdet' V2y 


detGi/2 Vdet'V2 


'E 

PL,PR 


p27rzTi (pop) — 2-7rr2 (p-p) 


(32) 


where c = 2, 


pop= pIGpl - PrGpr, 

P-P = pIGpl + PrGpr, (33) 

and 

Pl = ^ (G"^(n - Bw) + w) , 

PR=^{G~^{n-Bw)-w) . (34) 

The Laplacian determinant can be regularized as det^ = 
where the prime superscript indicates that the zero 
modes have been removed and p(t) is the Dedekind eta func¬ 
tion. 

Comparing with the standard formula for a CFT partition 
function on the torus one can read off the spectrum of scaling 
dimensions 


A = pIGpl + PrGpr + UrN^ + Y ’riR^n, 

nL>0 nR>0 

(35) 


where the last two terms correspond to the determinant of the 
Laplacian and represent, in canonical quantization, harmonic 
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oscillators indexed by positive integers nL,R- is the 

occupation number of oscillator 

For the SU(3)j^ WZW theory, we take G and B to be. 


1 

'2 

1 ' 

„ 1 

■ 0 l ' 

2 

1 

2 

’ ^=2 

-1 0 


Here G is proportional to the inverse of the Cartan matrix of 
SU(3). The above choice of parameters at the SU(3) point is 
consistent with that used in the numerical sections below. The 
relationship to the fields used in the previous section (note 
they were tilded) is simply a change of basis that diagonalizes 
G and rescales the diagonal elements to 1 i.e 

= h = \ ^{Xi+X2 ). (37) 

vStt V ott 

Since the terms in Eq. (21) correspond to the G term in the 
Hamiltonian (24) we deduce that the continuum version of 
the transformation i7su(3) ^ -ffsu( 3 ) + <lxQx is 

3 

GsU{3) Gsu{3) + <Ix^Gsu{3)- (38) 

This prediction will be tested with the help of accurate nu¬ 
merical calculations, discussed at length in the next section. 
Fig. 3 depicts the portion of the space that we traverse start¬ 
ing with our choice of parameters for the SU(3) model and 
varying G as a function of q^. The form of the G matrix in 
Eq. (38) is consistent with and expected from the Z 3 sym¬ 
metry, i.e., the conservation of Ni — Nq and N 2 — Xq - see 
Sec. II A. The Z 3 symmetry can be thought of as a ^ ro¬ 
tation in the root space of SU(3). In the effective field the¬ 
ory this is represented by the transformation on the currents 
(J“% J“^, J“^) (with superscripts defined in Eq. (18) ) as 
(J“i, J“^ J“ 3 ) ^ (J-“^ J“3, J-“i). In terms of the bo¬ 
son fields (j)i and (j >2 which live on the SU(3) root lattice, this 
amounts to 

In the X basis the Z 3 symmetry is represented by 

X^M' -X, M' (40) 

The G matrix in Eq. (38) is left invariant under the Z3 trans¬ 
formation. i.e. we have M'^GM' = G. One can show 
generally that any symmetric matrix left invariant by M' is 
proportional to Gsu{ 3 )- 

III. NUMERICAL RESULTS ESTABLISHING 
CORRESPONDENCE OF SPIN CHAINS TO CONFORMAL 
FIELD THEORY 

Having described the field theory for c = 2 spin chains, we 
now provide numerical evidence for the proposed correspon¬ 
dence. Our results first focus on various ways of calculating 


Im p = Im ^ 



Figure 2. Theory space of the two-component Tomonaga-Luttinger 
liquid with p = ^. The shaded region represents the “fundamental 
domain”; Because of the duality, different points in the theory space 
that are related by the duality are isospectral. The fundamental region 
is a set of representatives for all points related by the duality. I.e., 
starting from points in the shaded region, by mapping these points by 
the duality group, the entire theory space is covered. Some spectial 
points in the theory space are also marked: “SU(3)” represents the 
SU(3) WZW theory, and “SU(2)^” consists of two copies of SU(2) 
WZW theories, which may be realized, e.g, as two copies of the XXX 
spin chain. 


Im p 



Figure 3. Theory space of the two-component Tomonaga-Luttinger 
liquid with ^ = 1/2-1- \/2l2i. The point "‘qx = 0(SU(3))” corre¬ 
sponds to the SU(3) symmetric Lai-Sutherland model. The red line 
with arrows represents the points in theory space traversed as qx is 
increased from 0. (The arrows here do not indicate the renormaliza¬ 
tion group flow.) 

scaling dimensions (35), after which we discuss the proce¬ 
dure for extraction of the TLL parameters Gab and Bab for 
a,b — 1,2. We numerically confirm an important prediction 
of the field theory, namely Eq. (38). 

We carried out ED and DMRG calculations for periodic 
chains; finite size scaling of the energy gaps provides es¬ 
timates of the lowest scaling dimensions. For bigger open 
chains, we calculate the same information from the mutual 
information for spatially disjoint blocks. The mutual informa¬ 
tion measure is completely determined from the ground state 
wavefunction, making it useful for situations where obtaining 
excited states is difficult. 

The numerical calculations in this section were performed 
with a combination of our own codes and the Algorithms and 
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Libraries for Physics Simulations libraries [53]. 


A. Inferences from Exact Diagonalization and Density Matrix 
Renormalization Group 

For a ID periodic chain of length L, the scaling dimensions 
Xj, corresponding to the excited state with energy Ej, are 
given by, 


27rua;j 

— hif) — 


L L log L 


(41) 


where a is a model specific constant, v is the TLL velocity 
obtained from the finite size scaling of the ground state energy 
Eo, 


^ = e + _ - _ (42) 

L °° 6L2 + L2(iogi)3 1 1 

where e^o is the energy per site in the thermodynamic limit 
and c is the central charge and 6 is a constant. The form of 
the finite size corrections was derived by Itoi and Kato [40] 
for the SU(3) symmetric point (i.e. = 0); here we have 

assumed the same form holds for > 0- 

We note that the above formulas assume all excita¬ 
tions propagate with the same velocity v, while for multi- 
component TLLs, more than one velocity may appear in 
general. (For more generic models, these formulae need 
modifications; for example see the work of Ref. [54] on a 
SU(2) X SU(2) model.) In our model, a naive continuum limit 
and the bosonization analysis, (11) and (21), suggests that the 
excitations of the system, even when ^ 0, should be de¬ 
scribed by a single velocity. We will take this as our work¬ 
ing hypothesis. While our spectral analysis by ED/DMRG 
depends on this assumption, our later analyses based on the 
entanglement entropy and the mutual information do not. 



Figure 4. Velocity of the coupled TLLs as a function of q^. The 
inset shows the energy per unit length in the thermodynamic limit as 
a function of q^. The lines are guides to the eye. 

Figure 4 shows the TLL velocity v and ground state energy 
per site in the thermodynamic limit Coo as a function of q^ 


obtained by fitting our data to Eq. (42). Our results for the 
SU(3) symmetric point are in excellent agreement with ana¬ 
lytic results [55] and previous numerical studies [55-58]; for 
example, we get Coo = 0.29679 and v = 2.107(1) which 
are close to the exact results of 2 — In 3 — and 27r/3 re¬ 
spectively. Care must be taken in comparing our results with 
studies which parameterize the bilinear and biquadratic terms 
in the Hamiltonian 2a to be J cos 0 and J sin 9 with 0 = ^, 
thus requiring an additional factor of 1/y/2. We have used the 
value of the central charge c = 2, which we established inde¬ 
pendently from the scaling of the entanglement entropy (EE), 
discussed next. 

Before we proceed, we mention an important subtlety as¬ 
sociated with the choice of system sizes used in finite-size 
scaling. In a previous DMRG study on the SU(3) symmet¬ 
ric model. Ref. [57] showed the absence of the singlet ground 
state (scaling dimension 0 in the CET) for chains with lengths 
6M -F 2 and 6M 4-4, where M is a positive integer. [59] Thus, 
we restrict ourselves to analyzing chains with lengths that are 
multiples of 6. 


Central charge 

We establish the relevant region in parameter space where 
the TLL physics is expected to hold. Eor this purpose, we 
extract the central charge c, obtained from the scaling of the 
EE of a subsystem or "block", readily available in DMRG, as 
a function of its size 1. Eor open chains, the analytic form for 
the EE, denoted by S{1), is, 

5'(0 = -FS'o, (43) 

where Sq is a subleading correction. In Eig. 5 we show the 
profile of the EE and verify that the c = 2 fit to it is accurate 
for all qx > 0. [60] However, the EE profile has local structure 
occurring on the scale of three sites, that arise due to open 
boundaries. These are not captured by the leading term in 
Eq. (43). Other similar quality fits are possible with a lower 
value of c; we estimate c = 1.96 ± 0.05. Also note that Sq 
is non-universal; in this case dependent on q^ alone. This 
explains why the various curves in Eig. 5 differ despite having 
the same central charge. 

We pursue an understanding of the TLL behavior for all 
qx > 0 by considering the case qx ^ oo [37]. In this limit, 
the model is a purely classical one, with a macroscopically 
large number of ground states. To see this, we write out the 
Qx term on a bond {i, j) in terms of Sz and operators, 

Ul^Wx+h.c. = 2-3Sf-3Sf + ^sfsf + ^SlSi. (44) 

This expression indicates that the configurations 1 4 = 1± 1) and 
the configurations |±1 0) (and |0 ±1)) are exactly degenerate 
and have the lowest energies possible. This means that starting 
from a spin-1 "Neel" state, for example |4-1 — 14-1 — 1), one 
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Figure 5. Entanglement entropy as a function of block size (Z) and its 
fit to the formula (43) for c = 2 for a 150 site open chain for various 

(ix‘ 


can locally replace each | + 1 — 1) "dimer" by a | + 1 0) with¬ 
out changing the total energy. Thus, there is an exponentially 
large number of degenerate states. Adding the SU(3) sym¬ 
metric term lifts this degeneracy, but the model stays critical. 
Such a macroscopic degeneracy does not exist in the spin-1/2 
XXZ model in the Ising limit; this is why there is a finite value 
of anisotropy at which the spin-1/2 XXZ model ceases to be 
critical. 


Scaling dimensions and degeneracies 

In order to obtain multiple excited states in the same sym¬ 
metry sector (here sectors of definite Sz), we perform a state 
averaging procedure with two target states in the finite sys¬ 
tem DMRG method. A sequence of bond dimensions varying 
from m = 400 to m = 2000 states and periodic chains of 
lengths varying from 24 to 66 sites, were studied. For the ED 
calculations (from 6 to 18 sites), multiple excited states were 
calculated to give us a picture of the low energy degeneracy 
structure of this model. 

A note about boundary conditions is now in order. Work¬ 
ing with open boundary conditions, favorable for DMRG, can 
complicate the mapping of a spin chain to a conformal field 
theory: the notion of strict "conformal invariance" is broken. 
Hence we do not rely on open boundary conditions to give us 
a picture of the degeneracy structure of this system. That said, 
scaling dimensions can still be reliably numerically estimated 
from open chains. 

For the SU(3) symmetric model, it is analytically known 
that the first excited state is 18-fold degenerate in the "con¬ 
formal limit" and the second excited state is 16-fold degener¬ 
ate. However, in finite size simulations, the conformal limit 
is reached rather slowly as a function of system size; more 
specifically the lattice model (5) flows into the SU(3)j^ WZW 
critical point only logarithmically fast. Thus, we rely only on 
trends seen in the ED results. 


For the 12 site ED results, we observe that the low energy 
manifold consists of a non-degenerate singlet state, two sets of 
16-fold degenerate states [the occurrence of 16 being a con¬ 
sequence of SU(3) symmetry], followed by two degenerate 
singlets. As can be seen in Fig. 6(a), (d) on going from 12 to 
18 sites, the two singlets descend below the second manifold 
of 16 states: it is thus conceivable (though not rigorous), that 
these two states will join the 16-fold degenerate first excited 
states resulting in a 18-fold degeneracy in the conformal limit. 

Next, consider the effect of adding the term with > 
0. From ED, we find that the (exact) 16-fold degeneracy of 
the first excited state splits; the first excited state is now 4-fold 
degenerate, all corresponding to 5^ = 0 states, and the next 
excited state is 12-fold degenerate, corresponding to four sets 
of S'z = ±1 and two sets of Sz = ±2 states. Here too, the 
two degenerate singlets in the low energy spectrum descend 
to lower values on increasing the length of the chain, as can 
be seen in Fig. 6(b), (e) and (c),(f). Based on our experience 
with the SU(3) point, we conjecture that restoration of con¬ 
formal symmetry will result in the 4-fold degeneracy being 
transformed to a 6-fold degeneracy; although other possibili¬ 
ties are not completely ruled out based on this data alone. We 
expect this degeneracy structure to hold on varying only 
as long as the second excited state does not become the third 
excited state. 

The field theoretic prediction (38) confirms these infer¬ 
ences. Once the second scaling dimension exceeds the value 
of 1, which occurs around q^ « 0.5, there is a reorga¬ 
nization of energy degeneracies. For q^ ^ 0.5, we de¬ 
duce that the quantum numbers {ni,n 2 ,Wi,W 2 } [see Eq. 34] 
corresponding to the lowest 6 states are {±1, ±1,0,0}, 
{±1,0,0,0} , jo, ±1,0,0} and those for the next 12 
states are {^1, ±1, ±1, ±1}, {±1,0,0,±1} , {0,0, ±1,0}, 

{0,0, 0, ±1}, {0, ±1, Tl, 0} , {0, 0, ±1, Tl}. 

Figure 7 shows fits to Eq. (41), after taking logarithms of 
both sides, to extract the second scaling dimension X 2 , for 
various q^', similar trends are seen for the first scaling dimen¬ 
sion as well. The corrections to scaling are found to increase 
on going from the q^ — 0.0 to = 0.5. Whether these ef¬ 
fects are genuine deviations from the TEE physics or a lack 
of sufficient size to see "true scaling" can not be definitively 
established within our present methodology. We believe the 
deviations close to q^ « 0.5 are due to "energy crossings" 
(i.e. changing multiplet structure), causing additional level 
repulsions. Thus, one may need very large sizes to get precise 
estimates in this region. 

Despite this source of inaccuracy, the scaling dimensions 
vary within 10% when they are computed using Eq. (41) for 
fixed L, over the range of lengths considered (24 — 66 sites). 
The obtained values validate the correspondence between the 
lattice model and the GET and the general trends of their vari¬ 
ations with qx support our main conclusions. 
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Figure 6. Low energy manifold of the Hamiltonian (2a) for L = 12 (upper panels) and L = 18 site periodic chains (lower panels) for different 
values of q^. For = 0, shown in panels (a),(d) the first excited state is known to be 18-fold degenerate in the conformal limit. On increasing 
L, the inferred trend is that the two higher lying singlets (marked by red circles) descend to possibly join the 16-fold exactly degenerate states. 
(b),(e) and (c),(f) show similar trends for q^ = 0.2 and q^ = 0.5; in these cases the degeneracy structure in the conformal limit is narrowed 
down to a few possibilities. 



Figure 7. Finite size scaling of the second excited state energy gap 
for various q^;. Two independent DMRG calculations are performed, 
one each for the ground state in the Sz = 0 and Sz = 1 sectors. The 
fits to Eq. (41), along with the knowledge of the TLL velocity give 
the second lowest scaling dimension. 


B. Extracting the lowest scaling dimensions from mutual 
information 

It is difficult to target multiple excited states in DMRG for 
long chains, especially for a critical system where the entan¬ 
glement entropy grows logarithmically with system size. Thus 
it is extremely desirable to have a method to obtain scaling di¬ 
mensions that involves only the ground state. 

Typically this is achieved by measuring ground state corre¬ 
lation functions between two distant regions. However, in the 
most general setting, we a priori do not know the scaling op¬ 
erators on the lattice i.e. the operators whose expectations are 
to be measured. To obtain generalized correlation functions 
between two regions (say A and B) we calculate their com¬ 
bined reduced density matrix, for varying separations, and ex¬ 
tract the "mutual information" denoted by Iab and formally 
defined as, 

Iab = Sa + Sb — Saub, (45) 

where Sa,Sb, Saub is the EE of regions A, B and the union 
of A and B respectively. A schematic of the geometry used 
for this computation is shown in Eig. 8(a). 

The mutual information, unlike the block entanglement en¬ 
tropy, is not directly available in DMRG and must be cal- 
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A ^-► B 

jig.[ * □ 

(a) 



(b) 


For a finite periodic chain, one replaces Xij by the cord dis¬ 
tance L/tt sm{Trxij/L), this results in, 

rCC sm{Trx 32 /L) sin{TTX 4 i/L) \ 

3 \sin(7ra;3i/L) sin(7rx42/i)/ ' 

It is thus convenient to define the conformal ratio x as, 

^ ^ sin(7ra;32/L) sin(7ra:4i/L) 
sin(7ra;3i/ L) sin(7ra:42/ L) 

The notion of mutual information can be generalized beyond 
the von-Neumann entropy, which is assigned an index n = 1, 

fn) 

and thus denoted more generally by I\g. This is achieved by 
the following replacements in the CC formulae, 

1 “I” 71 

Si —> Sn, c —)■ —-—c. (50) 

bn 

Ref. [15] found that the true mutual information and the CC 
mutual information differ by a function f^{x), 


Figure 8. (Color online): (a) Geometry used to compute the mutual 
information consists of two 2-site blocks separated by distance r. 
The calculations were performed for a 150-site chain with r < 30, 
larger r data was discarded to avoid edge effects. Panel (b) shows 
F^{x) — 1, derived from the mutual information measure Eq. (45) 
as a function of the conformal ratio x (see Eq. (49)) calculated with 
DMRG. It was fitted to the analytic form, Eq. (53) to obtain the low¬ 
est scaling dimension. 


culated in a matrix product state (MPS) framework. (Practi¬ 
cally, this is achieved by reshaping all left and right optimized 
transformation matrices at the end of the DMRG calculation 
to get the MPS. Then, the reduced density matrix of disjoint 
regions is calculated using a partial-contraction scheme dis¬ 
cussed in Ref. [11]. More details of our calculations will be 
provided elsewhere.) The mutual information can also be cal¬ 
culated with Monte Carlo methods in sign-problem free sys¬ 
tems [13, 61]. 

We now discuss extraction of the lowest scaling dimension 
from Iab, for which we briefly present known results from 
the literature. To do so, we closely follow Ref. [15], whose 
notations we also use here. 

For a CFT, Calabrese and Cardy (CC) [1] argued that the 
entanglement entropy of two intervals A = [xi, X 2 ] and B — 
[x 3 , X 4 ] in an infinite lattice is given by, 

Saub = o log - + 2si, (46) 

3 V X 31 X 42 ) 

where x^ = Xi — Xj. The constant 2si is determined by de¬ 
manding that Saub — > Sa + Sb in the limit X 2 i,X 43 << 
2 ^ 31 , X 42 . Rewriting this formula in terms of the mutual in¬ 
formation (i.e. on subtracting out the single interval contribu¬ 
tions), one gets. 


tCC 

^AB 


3 V 2^31 2 : 42 / 


(47) 


I^AB-I^^^TB = f^^\^), ( 51 ) 

which is reparameterized as, 

= (52) 

n — 1 

Calabrese and co-workers [62, 63] have shown that for n > 
1 and in the limit of small x, 

- 1 = (£^)“«(») + (3^)'%4(») 

-F (higher order), (53) 

where a is twice the lowest scaling dimension Xi. The co¬ 
efficients S 2 (n) and 34 ( 71 ) are the contributions in the small 
X expansion coming from the two and four-point functions of 
the operator in the CFT with the lowest scaling dimension. 

Two concerns when using equation (53) in numerical sim¬ 
ulations are ( 1 ) it holds only for an infinite lattice and ( 2 ) it 
assumes that the non-zero contributions are solely from the 
operator with the lowest scaling dimension. However, for a 
finite system there are contributions from all operators. Thus 
the lowest scaling dimension fitted is simply an effective one 
trying to mimic the action of a linear combination of many 
(different scaling) operators. Empirically, for an open chain 
of 150 sites, all the errors (systematic and due to fitting) ap¬ 
pear to be within 10 %, which is roughly the error we also 
obtain from fitting to energies. 

Our results for fits to a power law for (x) — 1 for vari¬ 
ous qx are shown in Fig. 8 (b). The overall fits are reasonable, 
though there are local features not captured by Eq. (53): just 
like the case of the EE, these are attributed to open bound¬ 
aries. Such features are also seen in the spin-1/2 XXZ model, 
studied independently by Ref. [14]. 
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Figure 9. The main panel shows the first two scaling dimensions, x\ 
and * 2 , as a function of q^, obtained from finite size scaling of en¬ 
ergy gaps obtained from a combination of exact diagonalization (ED) 
and the density matrix renormalization group (DMRG) for periodic 
chains denoted by PBC. The c = 2 conformal field theory (CFT) 
prediction is also shown. The inset shows the lowest scaling dimen¬ 
sion obtained from the mutual information (MI) measure (see text), 
computed within the DMRG/Matrix product state framework for an 
open chain of 150 sites. 

C. Extraction of TLL parameters 

Fig. 9 shows the lowest two scaling dimensions obtained 
from finite size scaling of energy gaps as a function of q^. 
The inset shows the lowest scaling dimension from the mutual 
information method; with this metric, we were able to explore 
a larger range of q^. The general agreement (within errors) 
between these independent metrics confirms our that we can 
reliably calculate lowest scaling dimensions. Thus we proceed 
to discuss the extraction method for the four TLL parameters. 

Given a trial set of Gab and Bab, we calculate the lowest 
18 scaling dimensions, which need not be distinct, and denote 
them by Xj' . We then evaluate a cost function, 

(7(0,6) = (54) 

i 

and minimize it with respect to Gii,Gi 2 , G 22 and B 12 to ob¬ 
tain the best fit. We used the Nelder-Mead simplex algorithm 
built into the GNU Scientific library for this purpose. 

In order to confirm our inferences about the nature of the 
degeneracies in the low energy manifold, we attempted to 
fit to two degeneracy structures for the first and second ex¬ 
cited states. First, we assumed that the degeneracy (de¬ 
noted by gi) of the first two distinct scaling dimensions to be 
{ 91 , 92 ) = (6,12) and in the second case ( 51 , 52 ) = (4,12). 
In all cases, for q^ < 0.5, we found the former gave a signif¬ 
icantly better fit to the CFT formulae (35). In fact, attempts 
to use the (4,12) structure gave optimized solutions closer to 
a ( 51 , 52 , 53 ) = (4,2,12) degeneracy structure, hinting that 
the imposed structure was incorrect. The quality of our fits 
are checked by how well the scaling dimensions were repro¬ 
duced; for the correct degeneracy structure, these agreed to 


within ±0.03. 

The agreement of the values of the measured and expected 
scaling dimensions, shown in Fig. 9, strongly indicates an in¬ 
ternally consistent scenario for the lattice model to CFT map¬ 
ping. This is also equivalently seen in the extracted TLL pa¬ 
rameters, shown in Fig. 10, which are consistent with Eq. (38): 
they satisfy the expected relation Gu — G 22 = 2 ( 7 i 2 . The 
relative error in the scaling dimensions propagates to these 
parameters; for example, the overall error in Gu is roughly 
twice the error in X 2 - As expected from the duality explained 
in section IID, the scaling dimensions depend on B 12 up to 
an integer shift. Thus we focus on a particular representation 
and find that B 12 = 1/2 explains our data for all 5 ^,. Finally, 
even though we have shown data only for q^ < 0.5, the mu¬ 
tual information data in Fig. 9 (inset) suggests the validity of 
the theory for larger q ^. 



Figure 10. TLL parameters as a function of extracted from match¬ 
ing scaling dimensions to a c = 2 CFT. The main panel shows Gw 
as a function of q^ and we find Gw = G 22 = 2(7i2. The inset shows 
B \2 which is found to be constant. The CFT prediction is shown by 
the solid lines. 


IV. CONCLUSION 

In conclusion, we have developed an analytic correspon¬ 
dence between c = 2 free boson theories and microscopic 
spin-1 models, using bosonization techniques. For the par¬ 
ticular form of Hamiltonian considered (Eq. (2a)), we made 
a prediction for the value of the Tomonaga-Luttinger liquid 
(TLL) parameters as a function of q^, the parameter charac¬ 
terizing the lattice model. 

To build evidence on the numerical front, we performed ex¬ 
act diagonalization (ED) and density matrix renormalization 
group (DMRG) calculations to obtain the lowest scaling di¬ 
mensions from the energetics of the system; a scheme feasible 
for short periodic chains. However, our use of the mutual in¬ 
formation entropy between disjoint blocks, calculated solely 
from the ground state, provides a promising route to extend 
these calculations for long chains. 

























12 


Using this numerical data and the mapping from spin- 
chains to c = 2 theories, we deduced the value of all four 
TLL parameters as a function of q^. We expect our analy¬ 
ses to apply to more general situations, for eg. the model in 
Eq. (2a) with non zero qy. In future work, we aim to extend 
these ideas to calculate multiple low-lying scaling operators 
and dimensions of the c = 2 CFT using the correlation den¬ 
sity matrix [10, 11]. 

Our broader objective is an effort to develop generic meth¬ 
ods to map lattice models to multi-component field theories. 
We anticipate that this multi-scale modelling approach will 
be useful for understanding the physics at very large length 
scales; sizes that may not be directly accessible in numerical 
simulations. Once we have built confidence in the mapping 
between the lattice and continuum descriptions, we can use 
the (often known) predictions of the latter. 
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